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CH ■ Abstract 

We propose a computationally simple model of Auger and APECS line shapes from 
open-band solids. Part of the intensity comes from the decay of unscreened core- 
holes and is obtained by the two-body Green's function G-, as in the case of filled 
ON . bands. The rest of the intensity arises from screened core-holes and is derived using 

a variational description of the relaxed ground state; this involves the two-holes- 
one-electron propagator G U1 which also contains one- hole contributions. For many 
transition metals, the two-hole Green's function G w can be well described by the 
Ladder Approximation, but the three-body Green's function poses serious further 
problems. To calculate G u , treating electrons and holes on equal footing, we propose 
a practical approach to sum the series to all orders. We achieve that by formally 
rewriting the problem in terms of a fictitious three-body interaction. Our method 
grants non-negative densities of states, explains the apparent negative- U behavior 
of the spectra of early transition metals and interpolates well between weak and 
strong coupling, as we demonstrate by test model calculations. 

PACS Numbers: 79.20. Fv; 79.60.-i 

Keywords: Electron impact: Auger emission; Photoemission and photoelectron spectra 



1 Introduction 

The Core- Valence- Valence (CVV) Auger Spectra of solids with closed valence bands are well 
described [|l| in terms of the two-holes Green's function GL by the Cini [fj, |3| and Sawatzky Q 
(CS) model. 

The CS model allows to understand the phenomenology involving band-like, atomic-like and 
intermediate situations in terms of the U/W ratio of the on-site repulsion U to the band width 
W. For low U/W, the line shape is close to the self-convolution on the local one-hole density 
of states; with increasing U/W, the shape is distorted until, for a critical value of the ratio, 
two-hole resonances appear. They correspond to poles of Gu , as calculated from Anderson, 
Hubbard or related models. Atomic-like peaks and band-like structures are often observed ||T| 
in the same spectrum. For high U/W (atomic-like case) the Auger line shapes are so close to 
the free-atom spectra that they are labeled by LSJ terms and levels. Detailed studies of noble 
transition metals, like Au || and Ag || led to a very good agreement between theory and 
experiment, and also allowed the direct observation of off-site interaction effects 0. A similar 
success was achieved for covalently bonded solids ||. In all cases the diagrammatic expansion 

(2) 

of Gui is just a ladder of successive interactions between the holes without self-energy terms 
or vertex corrections. 

For open valence bands, the theory is much more complicated, and it is known from exper- 
iment that no such narrow atomic-like peaks exist any more. However for almost completely 
filled bands, when the number of holes per quantum state ii), C 1, remnants of the atomic 
multiplet structure are still observed and the closed-band theory can be extrapolated || . To a 
first approximation, one can assume that the valence electrons remain frozen during the core 
ionisation, and in the initial state of the Auger decay the valence configuration is the same as 
in the ground state. The intra-band shake-up effects are 0(rih) and can be accounted for by 



convolving with an asymmetric Doniach-Sunjic [10] line shape. In this way, the CVV spectra 

(2) 

are still described in terms of the two-hole Green's function Gw , computed in the absence of 
the core hole. Using Galitzkii's Low Density Approximation ITTJ (LDA) the dominant diagrams 
of the perturbation expansion of the G£, are just the same ladder diagrams which provide 
the exact solution for n^ — > 0. This gives a satisfactory explanation of the Auger spectra for 
nh « 0.1, which includes interesting cases like Ni and Pd [[HJ, even in finely dispersed form Jr3| . 
To extend the LDA to larger nh, a self-consistent version of the Low Density Approximation 
that involves using dressed propagators in the ladder series was naturally suggested J14J1 . How- 



ever an unexpected result cames from cluster studies |T5| which clearly demonstrated that the 

(2) 

ladder approximation with bare propagators is superior and is a good approximation to Gw 
for a wide nh range. This is the Bare Ladder Approximation (BLA). This approach has been 
useful to interpret the line shape of Graphite [ffHj. 



The next turn came from experiments on early 3d transition metals, like Ti and Sc |17j , 
that could not be interpreted by the above theory. The maximum of the line shape was shifted 
by the interaction to lower binding energy, which is the contrary of what happens in closed 
band materials. Qualitatively the CS model could work if one admitted that U < 0, and such 
an explanation has actually been proposed fll8| . 

However, no other evidence of U < was found; rather, it became clear that for almost 
empty bands one must formulate a new theory which is no simple extrapolation of the closed- 



band approach. Sarma |T9[ first suggested that the Auger line shape of Ti looks like some 
linear combination of the one-electron density of states and its convolution. Using this hint, 
and a general formulation of the Auger decay by Gunnarsson and Schonhammer |[20| , a simple 



explanation of the apparent negative- U behavior was found fzl|| . In this theory, the Auger 
line shape has two main contributions, that we call unrelaxed and relaxed, respectively. The 
unrelaxed contribution is obtained assuming that the Auger decay occurs while the conduction 
electrons are in their ground state \ip) in the absence of the core hole. The density of states 

(2) 

which shows up in this contribution is obtained by G w , like in closed-band systems. However, 
this is only a part of the story, and in many ways the easy part, because the screening electronic 
cloud that surrounds the core- hole can partecipate in the Auger decay. The ground state \(f>) in 
the presence of the core hole also enters the description. The relaxed contribution is computed 
with |0) as the initial state of the Auger decay. We summarise in the Appendix |A] the argument 
wich leads us to express the relaxed contribution in terms of a three-body (the Auger holes and 
the screening electron) density of states D; this is obtained from the Fourier transform of the 
t > part of the Green's function 

G(ai(3a h i l (3[a' l ;t) = (-z) 3 (^|T |4 ; (t) a^ (t) a 7i (t) a^a^a^ J \ip). (1.1) 

Here, operators are in the Heisenberg picture and we introduce a special notation ai,j3i... to 
indicate the set of quantum numbers of the local valence spin-orbitals belonging to the atom 
where the Auger decay occurs (the Auger site, for short). Experimentally, one can single 
out the relaxed contribution of ( |1 . 1| ) by properly fixing the photoelectron energy in an Auger- 



Photoelectron Coincidence Spectroscopy (APECS) experiment, |22|, |23|, ^4], |^| where the Auger 
electron is detected in coincidence with the photoelectron reponsible of the core hole creation. 
Fixing the photoelectron energy, the Auger electron measured in coincidence comes from the 
decay of a few dominant intermediate states in the presence of the core hole [2o| . 

Since ( |1.1| ) is hard to calculate and even much harder than of Gi, , we propose a simple 
approach in the spirit of the BLA. If we wish to design an affordable scheme in this difficult 
problem, we must be prepared to adopt a series of approximations: Cini model, Ladder Ap- 
proximation and a new one, which we call the Core Approximation. Section [^ is devoted to 
the fomulation of the new scheme. The degree of validity of our approach will be investigated 
by comparison with exact results from cluster calculations in Section [3[ 

2 Perturbation expansion of the two and three-body 
Green's functions 

In the Cini model, we express the Auger line shape in terms of local Green's functions, like the 
one of Equation (|1.1|); moreover, we calculate them by considering only the local scattering at 



the Auger site, since in this way we drastically simplify the algebra, and the line shape is little 
influenced by scattering at the other sites |^7f| . Many Auger line shape calculations in solids 
have been performed in this way. 



The local interaction approach has been extended successfully to open bands [12|, |13| . In 



the present work we want to study its application to more general n^. This should allow to 



extend the analysis to several transition metals, giving at least a qualitative understanding of 
their spectra, which is currently a difficult task. 

In order to properly evaluate the results one should bear in mind that currently even for 

(2) 

Gu we have reliable recipes only for n^ less than « 0.25. This problem involves one more 
body and highly excited states of interacting systems; consequently even the main features of 
the solution are quite an unsettled question. 

2.1 Ladder Approximation to the two-hole propagator 

(2) 

The time-dependent Green's function G t allows the usual perturbation expansion, in terms 
of time-ordered products of interaction-representation operators, namely: 

G( 2 >(a,A,#ai;*) = X)H) n / <**!••• / dt n (T[ai i {t)a\ i {t)H u {t l )...H u {t n )a p[ aA) c , 

n J —oo J— oo 

(2.1) 

with the average taken over the non interacting ground state \ipo) and the sum restricted to 
the topologically inequivalent, connected diagrams. Ordering the spin-orbitals in an arbitrary 
way, we may write the Coloumb valence- valence interaction between those of the Auger site in 
the form 

Hu= Yl U wiPin a tati a n a Pr ( 2 - 2 ) 

m< v upi<n 



The diagrammatic method develops Equation (|2.1| ) in terms of local non-interacting time or- 
dered one-body propagator So (atj, Ai t) 

S (a h A; t) = S% K A; t) - S e (A, a,; -t) , (2.3) 

where 

S%(a h [3 l ;t) = -ie(t)(at tl (t)a Pl ), (2.4) 

S e (A, a,; -t) = -id (-£) (a A < (t)>. (2.5) 

here the average is taken over the non-interacting ground state \ipo) with energy Eq. For 
n h < 0.25 and a range of U/W, a good start is provided by the BLA. In the BLA, one selects 
the series of ladder diagrams which are free of self-energy insertions and vertex corrections. 
This is equivalent to the approximate factorisation 



(T {< (t) 4, (*) Hu (h) ...Hu (t n ) a p[ a a ^) c « ]T U, imPlTl 

tH<Ul,pi<Tl 

[(T {< (£) al (t) a n {t x ) a Pl (*0})<T {< (h) a\ (h) H v (t 2 ) ...H v (t n ) a%a a[ }) c ] . (2.6) 
Since 

g (ai/3 b Tipi;t - t x ) = (Wo*, it) a\ t {t)a n (t 1 )a Pl (ti) }>, (2.7) 



is the non-interacting Green's function, the BLA leads to 
G {2) (aiPi, /3,'aj; t) = g (a/A, M; t) 



^ U^ VlPlTl J dtxg(ai(3 h Tipi\t-tx)G { - i) (iJLiPiiPla'iiti). (2.8) 
H<n,pi<n 



Going to frequency space, this becomes a linear algebraic system. Equation (|2.8|) is called Bare 
Ladder Approximation because it uses undressed single particle propagators. It is the exact 
solution for n^ = 0, and remains a good approximation for a useful range of n^. It is simply 
equivalent to keeping only those diagrams that remain in the closed-band limit. For closed 
bands, the one-body propagators of Equation (|2.3|) reduce to the S$ part; therefore, lines that 
start as hole lines never go back in time and remain hole lines thoughout the diagrams. So, the 
order of times t, ti, t 2 , ■ ■ ■ ,0 remains fixed (decreasing) despite the presence of the time-ordering 
operator T. 

The convolution form of Equation ( |2.8| ) has furter important consequences. Mathematically, 
it is a Dyson equation in which the U matrix is an instantaneous self-energy Therefore, it grants 
the Herglotz property: for any interaction strenght, G^ generates non- negative densities of 
states. The Herglotz property is a basic requirement for a sensible approximation, yet it is not 
easily obtained by diagrammatic approaches. Its achievement is one of the most interesting 
features of this simple approximation. 

2.2 Ladder Approximation to the three-body propagator 

The Green's function (|1 . 1|) yields the expansion 

/oo /»oo 

dtv dt n 

-oo J —oo 

(T {4 ; (t) a\ (t) a 7i fa) H v fa) ...H v fa) a\ { a p[ a a ^) c . (2.9) 

This describes the propagation of two-holes and one electron in the final state, or, if the electron 
and one hole annihilate, a one-body propagation results (see Section |2.5| ). In proposing an 
approximation to Equation (|2.9|) , we may proceed by analogy with G 2 U . For G^ it is natural to 
propose a bare ladder approximation in which the lines that start as electron (hole) lines never go 
back in time and remain electron (hole) lines throughout the diagrams; such an approximation 
should have essentialy the same physical contents as in the previous case and work properly in 
a wide range of U and fillings. However, the series cannot be summed easily like in Equation 
Q2.8|) , because with three bodies involved we meet an extra difficulty For instance, let Hu (ti) 
produce an interaction between the two holes in a given term of the expansion; then, the 
electron line overtakes time ti] therefore the diagram does not yield a convolution of a function 
of t — ti times a function of ti. This undesirable feature can be removed by using the identities 

S%(a h p l ;t) = tJ2So(ai,r,t-t')S%( 7 ,(3 l ;t>). (2.10) 
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Figure 1: Second order contribution to the three-body Green's function. Using Equations ( 2.10J ) we cast it 
in the form of a product of three "blocks" . These are easly dealt with by a Fourier transform. 



and 



S e fa, A; t) = i J2 So K r, t - S e (7, A; • 



m 



where the summations run over all the complete set of spin-orbitals. These identities are derived 
in Appendix |B[ Here we note that in the limit t = t' — ► + , we get 



So(ai,f3i;t) -> -iy^^ala-f) (a\ap) , 



(2.11) 



and 



So {an, Pi; t) -*i) j (a ai al)(a 7 al i ). 



(EH) 



Since the ground state is not the hole vacuum, (a 7 ajj) 7^ S(j,j3); so the correct limits are 
obtained only thanks to the completeness of the 7 set. 

To see the use of Equations ( |2.10|) , consider for instance the application to one of the 
second-order contributions to ( |2.9| ). Using the standard diagrammatic rules, the l.h.s. of the 
uppermost Equation in Figure (|l|) reads (understanding a sum over repeated indices): 



G 2 (a l [3 nh7 ' l P' l a' l ;t) = -U l 



mnpin^ p-'i v [p'it[ v % ) 




dt\ dti 



[S e (7,, v x ; 1 - ti) S% (A, Pf, t - h) S f l fa, p\\ t - t 2 ) 

S h fa, r{\ ti - h) S e (71, t(; h) S h Q {u[, A'; h) S% fa, a[; t 2 )] , (2.12) 



Using the identities ( |2.10| ), one casts Equation ( |2.12| ) in the convolution form 



P /*00 

^ / / dtx dt 2 [Go (aiPryi, Pipia; t - t ± ) G (a/im, 77/pJ; t\ - t 2 ) G (^v'fi, iAa'u h)] ■ 

(2.13) 

where 

G (a/? 7 , 7W; t) = S e (7, Y; *) S# ((3, (3>; t) S* (a, a'; t) . (2.14) 

This is the expression, that we represent pictorally as the r.h.s. of the uppermost Equation 
in Figure (|I|). In this way, introducing a fictitious x interaction vertex, along with the true 



interaction vertex (dot), we obtain a factorisation similar to Equation (|2.6| ) and the diagram is 
cast in the convolution form. This useful property extends to all the diagrams of the bare-ladder 
approximation. 

2.3 Core- Approximation 

We achieved Equation ( |2.13| ) in the form of a convolution, which is a big simplification. However, 
the infinite summations (one for each x interaction) are a high price to pay for that. They 
arise because Equations ( |2.10| ) imply a summation over the complete set of 7. On the other 
hand, since we use a local Hu, the dot interaction involves only local matrix elements between 
spin-orbitals; so, we are only interested in the local elements S h,e (ai,/3i;t). Physically, we may 
expect that only the sites which are closest to the Auger site give an important contribution to 
the summations, and we can actually work with a limited set to express the local S h,e (an, /3f, t) 
elements. Larger sets will lead to more precise results, at the cost of more computation. Here, 
we wish to explore the simplest approximation, by drastically limiting the 7 summation to the 
local states 7/. We have observed above that summing over the complete 7 set is necessary to 
get the correct zero-time limit. In turn, this is essential to preserve normalisation. If we simply 
replace 7 with 7; in Equations ( |2.1U| ) this condition is violated (see Equations ( |2.11| )); therefore 
we introduce a set of functions R (ai,(3i;t) which approximately factor the propagator in 
analogy with ( |2.10|) according to the ansatz 



H)« (*)«/%> « J2 R+ («i,7i;*~0 « (0a A ), (2.15) 

and 

(-i)(a ai (t)al) w J2 R ~ ^ hll ;t-f) (a 7i (*>£,), (^TBJ) 

where t' is any time intermediate between and t; the R^' functions are computed for any t 
by solving the system 

(-i)(a) ai (t) ai3l ) = ^R + (aulut) (a\a Pl ), (2.16) 

11 



and 

(-i)(o a; (t)o^) = ^R~ {a h -fi-t) {a^al). ( gTEf ) 



a 



The system must be identically satisfied for any t and in particular the correct t — > + limit is 
granted. 

Equations (|2.16|) are correct in the limit of core states, when S h ' e (oti, /3f,t) is diagonal in 



its indices and coincides with R ±m , therefore we call this the Core Approximation (CA). The 
ansatz is also correct in the strong coupling case, when localised two-hole resonances develop. 
This is appealing, since the strong coupling case is the hard one, while at weak coupling 
practically every reasonable approach yields similar results. Thus, we regard the ansatz ( |2.15|) 
as a physically motivated aproximation, which must be tested against exact results for its 
validation. 

2.4 Summing the Three-body ladder 

Working out the Core Approximation (CA) like in the example fl2.13|) one can compute all kinds 
of ladder diagrams, to all orders. The partial sum of the series ( f2.9| ) that one obtains in this 
way will be referred to as Core-Ladder- Approximation (CLA). From now on only local indices 
appear so we shall dispense ourselves from showing this explicitly. In the exact expansion (|2.9| ), 
with the local Hy denoting the interaction Hamiltonian ( |2.2| ) between fermions at the Auger 
site, the outgoing holes at time t are labeled a and (3 (creation operators) and the outgoing 
electron 7 (annihilation operator). In the CA for the n-th term of 



/oo /»oo y*oo 

dh / dt 2 - ■ dt n J2 ^W 

00 J' 00 J' 00 fJ,<V, P <T 

(T |4 (*) ajj (t) a 7 (i) oj, (t x ) 4 (i x ) a T (t x ) a p (t x ) H v (t 2 ) ...H v (t n ) a+^/cv j) c , (2.17) 

we consider the times ordered with t > t\ > . . . > 0; let us first analyse the contribution 
which arises when Hu(t\) describes hole-hole scattering. In the expansion of the T-Product of 
Equation (|2.17|) by Wick's theorem we contract the a and (3 creation operators of the outgoing 
holes with the annihilation operators at time t\, obtaining among other terms: 



(T {4 (t) a p (ti)})(T |oJ (t) a T (*i)})(T [a 7 (t) oj (t x ) 4 (t x ) H v (t 2 ) . . . a\,apa a , ]) c , (2.18) 

This, however cannot yet be factored in the convolution form because the last term of Equation 
( |2 . 1 8|) depends on t. In the ladder approximation, the annihilation operator a 7 (t) can be 



contracted with any creation operator at whatever time t 2 , ■ ■ ■ ,t n or 0. Using the CA we can 
insert a x vertex at time t x on the 7 electronic line (like in Figure ([I])); namely 



(T {4 (t) a p (£i)})(T |oJ (t) a T (£1) j)(T |o 7 (t) a\ {t x ) a\ (t x ) H v (t 2 ) . . . a\,a p ,a a ,j) c = 

(T{ai(t)a p (t 1 )})(T{al(t)a T (ti)})^ R ~^^^ t - t ^ 

t 

(T |a € (h) 4 (h) 4 (ti) H v (t 2 ) . . . a\,a p/ a a ,}) c . (2.19) 
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Figure 2: Diagrammatic representation of the contributions of order n to the CLA (Equation (2.24)). The 
fictitious x interaction vertex represents a R function as explained in the example (2.13). 



This result can be easily verified by expanding with Wick's theorem, and then applying Equa- 
tions fl2.10| ) although the proof is somewhat lengthy. In Equation ( |2.19| ), the propagation 
between times t and t\ is described by the factor 



G (a/37,£rp; t - h) = FT (7, £ t - h) S% (/3, r; t - h) S% (a, p- 1 - h) 



(2.20) 



this is similar to a non-interacting three-body propagator, except that the S e has been replaced 
by a R~\ as a shorthand notation we underline the electron indices that correspond to a Rr 
factor. In a similar way Hjj(ti) can give other possible ladder contributions to Equation (|2.17f) ; 
these come from the interaction of the 7 electron with one of the two holes; in this case is the 
other hole that goes to a x vertex and we meet the factors 



G (a/? 7 , V/3V; t - h) = S e (7, 7'; t - h) R+ ((3, /?'; t - h) S% (a, a'; t - h) 
Go (aj3% V/9V; t-t 1 ) = S e (7, 7'; t - h) S% ((3, /?'; t - tx) R + (a, a'; t-h) 



(2.21) 
(2.22) 



where underline hole indices correspond to a R + contraction. 

A glance to Equation ( 2.17Q reveals that the last factor in the rhs of Equation ( |2.19|) is the 
matrix element which enters G n _i; hence we obtain the CLA form for one of the term of order 
n of the development, corresponding to the diagram (a) of Figure (H) 



U„ 



/ j K I'I'I'T 

fl<V,p<T 



E 

L (, 



dtxGo (a/?7,^rp;t-ti) G n _i (/j,u£,'y'(3'a';t 1 ) 



(2.23) 



In other terms of the expansion, Huiti) describes the interaction of the 7 electron with 
the a and j3 holes; in addition, each contribution has an exchange counterpart, obtained by 
crossing the hole lines going from t to t±. Thus, applying the CA to the other contractions of 
eq.( 2.17 ) we get the six contributions represented in the fig.(Q). Letting n — > 00 and Fourier 



transforming we obtain a linear system of equations for the CLA form G CLA (a/37, 7' '(3'a!\ uj) of 
the three-body Green's function, namely 



G CLA (a/37, Y/?V; uj) = B (a/57, 7^'a'; w) - 

J^ U^prBo (a/37, §rp; uj) 




p<T 



J2 ^W l B o ("07, ^P! w) G CLA O^r, 7W; w) + 

A, (a/3 7 , pr£; w) G CLA (£/p, Y/3W; w)] } . (2.24) 



fl<V,p<T 



where the B n functions are 



Bq (a/37, l'P' a 'i u ) — Go (a/37, l'P >a> i u ) ~ Go ((3aj, j'fl'a'] uj) , 

#o (a/fy, i r P; w) = G o (a/97, £rp; uj) - G (a/3% £pr; uj) , (2.25) 

B (a/37, ^fp; w) = G (a/3_7, i/£p; w) - G (a/37, i/p£; w) , 
5 (a/37, pr£; w) = G (a/37, P^; w) - G (a/3_7, p£r; w) . 



The second, third and fourth lines of (|2.24 ) come respectively from the (a,b), (c,d) and (e,f) 



diagrams of Figure (|2|); while the first two contributions come from the hole- hole interaction 
the others come from electron-hole interactions and convey information on the screening effects 
due to the electronic cloud which forms as a response to the deep electron ionization. 

2.5 Single-particle contribution 

At the same level of approximation we must consider the case when one (or both) the holes 
produced by the Auger transitions has the same spin as the screening electron. Consider the 
spin-diagonal components G (a/37, 7' (3' c/;t), (with cr a > = a a and so on); when the hole /3 has 
the same z spin component as the electron, then, contracting at it) with a 7 (t) and a pi with a , 
in (|2.17|) one obtains the extra contribution 



G sp (a/? 7 , 7 / /3'a / ; t) = (-1) (o^a 7 ) (a\,api)S (a, a'; t) . (2.26) 

where S (a, a'; £) stands for the time-ordered dressed one-body Green's function that can be 
expanded with 

/oo ^00 

dtv dt n (T{ai(t)H u (t 1 )...H u (t n )a a i}) c . (2.27) 

-00 J —00 



and summed with the Dyson's equation [28| in terms of proper Self- Energy diagrams. 




Figure 3: Topologically distinct contributions to the Second order self-energy. No tadpole diagrams occur 
because we expand in terms of Hartree-Fock spin-orbitals. 



At second order of Q2.27 ) we find eigth diagrams free of Hartree-Fock (tadpole) insertions 
that give rise to four topologically distinct ones as shown in Figure (^), which shows how the 
second-order self-energy can be expressed in terms of the hole-hole-electron Green's function 
G. The rules are: (i) write down a diagram for G; (ii) join the incoming electron line with one 
of the hole lines into an interaction vertex and draw an interaction line to the other hole; (iii) 
repeat the operation with the outgoing lines; (iv) do all that in all possible ways. This is a 
general result that can be extended to all orders yielding the Dyson form for the S (a, a 1 ; uj) 



S (a, a'; uj) = So (a, a'; uj) 



E E 

p! <U' ,p' <t' fJb<V,p<7 



U fivprU ii'v' p'r' 



{[So (a, p; uj) G (pur, v'r' ' p'\ uj) — Sq (a, r; uj) G (pup, v'r ' p'\ uj)} S (p , a; uj) + 

[So (a, r; uj) G (pup, p't'p'; uj) — S (a, p; uj) G (pur, p'r'p'; to)] S (u 1 , a'; uj)} ; (2.28) 

at the second order this reduces to the four terms of Figure ([3]) . We get a conserving approx- 
imation for the proper self-energy, using the CLA approximation for the three-body Green's 
functions of Equation 



By solving Dyson's equation (|2.28|) one can model XPS spectra from valence bands with 
low band filling; this is another field of application of our approach. 

Having computed S, one finds how the Fermi energy Ep is renormalised by the interaction. 
In general, we shall reference all the three-body spectra to the value of Ep which are obtained by 
the corresponding one-body calculation. In the next Section, when comparing the line shapes 
with exact diagonalisation results, we shall align the Fermi levels accordingly. 

The deep hole attracts a screening electron that can be directly involved in the Auger decay; 
this is the physical origin of the contribution (|2.26|) to the three-body Green's function. Locally, 
such processes leave the system with one hole in the final state. The presence of an one-body 
contribution in the Auger spectra from transition metals like Ti or Sc has been pointed out 
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already [ |I9| , |21|j . Besides the three-body ( 2.24|) and one-body Q2.26| ) diagrams, there are mixed 



contributions. Physically they represent interference contributions in which the system evolves 
from one-hole states to two-hole one electron states and back. Here we are going to neglect 
those diagrams for the sake of simplicity. Such terms can only be significant if the interaction is 
strong enough; however for strong interaction the one and three-body contributions are widely 
separated in frequency, so we may argue that the interference terms are not very important in 
general. Thus, we write 

G {a(3j, i(5'a'\ u) = G CLA {a(3j, i(3'a'\ u) + G sp (a/3 7 , i$ot\ u) . (2.29) 

3 Comparison with exact diagonalisation results 

In this Section, we wish to test the CLA results against those of a model system that can be 
diagonalised exactly. 

3.1 Model Cluster 

We consider a 5 atom cluster, with 2 levels for each atom; the one-body basis elements are 
\sia) with s = 1, . . . , 5 the site index, i = 1, 2 the level index and a =|, j for the spin direction; 
the one-body energies are denoted by e S j. The atoms 1 . . .4 occupy the vertices of a square, 
and the Auger atom is at site 5 above the centre. Although this is just a model, we wish to 
use it to outline a possible procedure for the analysis of actual experimental data. Suppose we 
know a simple one-body Hamiltonian H Q of the system, of the same level of sophistication as a 
tight-binding model of a solid, with a nearest neighbor hopping term: 

H o = Yl Usi ^ si + Yl Yl T if' a ii« a s'i^ ( 3 - X ) 

sia ija (ss') 

with real Tf s . For nearest neighbor sites we assume identical matrices 

Next, suppose we have an estimate of the Hubbard U and we wish to model the interactions 
by 



Hffub — U 



} J n s iin S ji + 22 n sian s2 a 



*'./ 



(3.3) 



The full interacting model Hamiltonian cannot be taken to be the sum of Hq and Hn u b, because 
Hq must contain the effects of the interactions on the one-body states, roughly, at the Hartree- 
Fock level. A mean field average of the interaction must be subtracted out. Therefore, we 
assume the model Hamiltonian 

H = H Q + H l , (3.4) 
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where 



and 



Hi = H Hub — Vj 



H-Fi 



(3.5) 



V H -f = U < ^2 [(n S i-\)n sji + n si1 (n sji )} + ^ [( n sia)n s2a + n s i a (n s2a ) + 

(. sij sff 

ali (T a s2cj {a s i a a\ 2(T ) + h.c.j j. (3.6) 



where mean values are taken over the non-interacting ground state. In principle, if H were 
a Self-Consistent-Field Hamiltonian, a self-consistent calculation of e parameters and mean 
occupation numbers would be a more precise procedure. In this way, the Hartree-Fock con- 
tribution to the self-energy would be automatically enbodied in the bare propagators. In the 
present paper, we use non-interacting ground state averages, which is a cruder approximation, 
because we wish to privilege the simplicity of the procedure; the simplified approach already 
gives gratifying results (see below). 



3.2 Exact diagonalisation calculations 

We have chosen the following parameters in eV: e s i = and e s2 = 1 for all sites; tn = 1.4 and 
ti2 = 0.4. In this way the bandwidth W, that is the difference between the extreme eigenvalues 
of Hq, is ~ 9.6 eV, a reasonable value for a transition metal. We want the densities D which 
are related to the positive time part of the Green's functions ( P-.lfl , and we are particularly 
interested in the diagonal elements D (5z<Tj, 5jaj, hkok\ w); these are usually the largest and 
must be non-negative. 

These are given by (dropping the site index s — 5) 



D(ia u jaj,k(Tk;w) = J^ 






s 



(w + E%l x - E 



(3.7) 



where N is the total number of electrons in the cluster, 



is the interacting ground state 



with energy E; \E^_^) is the complete set of interacting eigenstates with (N — 1) electrons. In 
one set of calculations we start from a ground state with N = 18 valence electrons, since the 
maximum occupation in the cluster is 20, we have two holes. We assume a non-magnetic ground 
state, with one hole for each spin; the configurations are ( 2 ) = 100. In the final state, with 
two Auger holes of opposite spin and a spin-up electron we find ( 2 ) ( i ) = 450 intermediate 
states with (N — 1) electrons. 

In another set of calculations we start from a ground state with N = 16, that is, two holes 
for each spin, and get \ip) and E by diagonalising a ( 2 ) x ( 2 ) = 2025 x 2025 matrix, while 
for |.Ejv-i) we mee t the maximum size of matrices in our calculations, namely ( 2 ) ( 3 ) = 5400. 

The site s = 5 was chosen to be the Auger atom because in this way the problem is highly 
symmetric, and we can classify all the states according to the Irreducible Representations {IR) 
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of the Ci v Group, the Group of the square. Carring out in Equation (|3.7| ) the summation on 
the IR of the group we obtains 

D fajaj, ka k ; w) = £ E | (MJ^ka.^f) |* 5 (u + E^f - E) , (3.8) 

IR n 

but the operators a ia belong to the totalsymmetric representation A\ of and so \ip) and |-E]v_i) 
must belong to the same IR. Since the ground state \ip) for N = 1Q and iV = 18 is total- 
symmetric we could restrict the summation to the A\ sector. 

D fa, fa, ka k ; uj) = J2 | {MaA°^ k \ E N-?) f 6 (« + E< N-f - E ) ■ (3-9) 

n 

with a strong reduction of computing time. 

3.3 Comparison and evaluation of the CLA 

By the cluster calculations we want to study the physical consequences of the hole-electron 
interaction in the three-body density of states. Also, we wish to estimate the degree of validity 
of the CLA. We performed the comparison with the exact results for several values of U/W. 
By considering various possible electron populations N < 20 in the cluster we also varied the 
mean occupation of the Auger site in the non-interacting ground state 

(n) = - ^ ( n 51<x + n 52a) ■ (3.10) 

a 

The test becomes more severe when (n) is reduced towards half filling and U/W is increased. 
Many different densities of states are obtained from the matrix elements of G(u); in the context 
of our theory, the density 

D lh (uj) = (^l4iT a 52| a 5iT 5 (w + H-E) 4 1T a 52i a 5 i T |^). (3.11) 

is of special interest because in the diagrammatic series the annichilation of the spin-up electron 
by the hole of the same spin is particularly strong so we may expect that the one-body term is 
important. By contrast, in the density 

D 2 hie{u) = (V'l4iT a 5i| a 52T^ (u + H-E) 4 2T a 5U a 5 i T |V>). (3.12) 

the same contribution should be smaller and possibly absent. We shall consider these two 
examples in turn. In Appendix |Cj we detail the application of the CLA to the problem at hand. 
Although the analytic development can be somewhat boring, the maximum size of the matrices 
involved is just 8. 

The left frame of Figure (|4]) shows Dih for U/W=0.25 with N = 18, which yields a popula- 
tion (n) = 0.86 on the Auger site. The density is dominated by a single peak at binding energy 
~ 4eV, but also shows a pair of wings. Since correlation effects are moderate and the filling is 
fairly high, the CLA is still in good agreement with the exact results. 
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Figure 4: Binding energy dependence of D lh for (n) = 0.86 with U/W=0.25 (left frame) and U/W=l (right 
frame). Heavy line: exact result; light: CLA; dashed: 1-body contribution to the CLA result; dotted: 3-body 
contribution to the CLA result. The line shapes have been convolved with a Lorentzian (FWHM=0.75 eV). 



In the right frame of Figure (|]) the interaction is increased to U/W = 1. The pattern does 
not show any major changes, except a shift and an increase of the structure at high binding 
energy. The shift (~ 5 eV) is large, which is understandable because the filling is high and the 
screening uneffective. The value of U is well outside the scope of weak coupling approaches but 
the CLA still reproduces the exact results rather well. 

Both the line shapes and the accuracy of the CLA are sensitive to (n). In a series of 



calculations, we set N = 16, which yields (n) = 0.72. 





Binding Energy (eV) 



Binding Energy (eV) 



Figure 5: Binding energy dependence of D lh for (n) = 0.72 with U/W=0.25 (left frame) and U/W=l (right 
frame). Heavy line: exact result; light: CLA; dashed: 1-body contribution to the CLA result; dotted: 3-body 
contribution to the CLA result. The line shapes have been convolved with a Lorentzian (FWHM=0.75 eV). 
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Figure 6: Binding energy dependence oi D 2 hie for (n) = 0.86 with U/W=0.25 (left frame) and U/W=l (right 
frame). Heavy line: exact result; light: CLA; dashed: 1-body contribution to the CLA result; dotted: 3-body 
contribution to the CLA result. The line shapes have been convolved with a Lorentzian (FWHM=0.75 eV). 



In the left frame of Figure (pi) a moderate interaction U/W = 0.25 is assumed, and the 
exact results show a threshold peak centered at the Fermi level accompanied by a rich structure 
between 5 and 11 eV. The CLA underestimates somewhat the intensity of the main peak and 
does not faithfully reproduce the shape of the high binding energy structure, however the 
position of the main peak is correct and the overall line shape is in fair agreement. Moreover, 
the performance of the CLA does not break down quickly with increasing U as weak-coupling 
approaches tend to do, but remains fairly stable. This can be seen in the right frame of Figure 
(|JP, where the comparison is done with U/W—l. Due to the effective screening in this case, the 
increase of U does not cause a dramatic shift of the high binding-energy structure; this is borne 
out by the CLA and we may still speak of semiquantitative agreement. The CLA also explains 
the increase of the relative weight of the 3-body contribution with reducing band filling. 

As one could expect, D 2 hi e (u) has much more weight at high binding energies than D lh (u), 
as one can see in Figure (|6|); in the right frame we have chosen the moderate coupling case 
U/W=0.25 with (n) = 0.86, like in Figure (|). 

The overall agreement of the CLA with the exact results is quite good for the main peak, 
and even the small structures at lower binding energy are reproduced in some detail. With 
increasing U, these structures are largely lost (right frame of Figure (Q)); the CLA deteriorates 
in those regions where the intensity is low, but the overall line shape is still satisfactory. 

The most remarkable feature is that the large increase of the interaction strenght produces 
a rather small increase of the binding energy of the dominant peak compared to the left frame 
Figure (||). This would not be understandable if the peak were a two-hole resonance. 

Physically, the shift is refrained by the screening of the two holes by the electron, which also 
becomes more effective with increasing U. Our simple approximation correctly accounts for the 
screening effect and consequently predicts the position of the main peak quite accurately. It is 
not surprising that L>2hie(^) is more difficult to approximate in detail than Dih(uj), since many 
more, highly excited final states are reached with three quasiparticles than with one. 
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Figure 7: Binding energy dependence of D 2 hie for (n) = 0.72 with U/W=0.25 (left frame) and U/W=l 
(right frame). Heavy line: exact result; light: CLA. The line shapes have been convolved with a Lorentzian 
(FWHM=0.75 eV). 



A particularly hard case is shown in Figure fl7|) where the filling is rather low ((n) = 0.72). 
No single-particle contribuition exists in this case because (V'l a 5iT a52 Tl^') = 0- 

For U/W— 0.25 (left frame) the agreement is fairly good. For U/W=l (right frame) the 
exact line shape shows three partially resolved broad peaks covering a range of ~ 10 eV in 
binding energy. The CLA in such severe conditions misses the main peak position by ~ 1 
eV and underevaluates the low binding energy shoulder; however, we can still claim at least a 
qualitative agreement with the results of the exact calculation. 

In all cases we find that the Herglotz property is fully preserved; this is a most valuable 
feature which is not easily obtained for approximate three-body propagators. For instance, the 
approach of Ref. |H|] fails in this respect at strong coupling. 

Comparing the two frames of Figure (0) we observe the apparent negative- U behaviour: 
increasing the interaction U, the main peak shifts towards lower binding energies; this is a 
consequence of the interaction of the screening electron with the two Auger holes. We remark 
that a high enough n^ is necessary to build up a localised screening cloud. This is why the 
negative- U behavior is observed in the early transition metals, but not in the late ones. 



4 Conclusions 

To model Auger and APECS line shapes, we need a reliable practical recipe for a 3-body prop- 
agator, which poses much more difficult problems than the 1 and 2-particle Green's functions. 
The Core-Ladder-Approximation that we are proposing produces a dramatic simplification. 
The problem is reduced to the calculation of the mean-field one-body Green's function and to 
the diagonalisation of small matrices. Comparison with exact model results shows that one 
can capture the essential physics of a formidable problem which otherwise would require an 
excessive computational effort. 
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Experience with the 1 and 2-body problems suggests that the 3-body ladder series, without 
self-energy and vertex corrections, is a prosiming approximation: it treats electrons and holes in 
the same way, and produces non-negative densities of states. From a mathematical viewpoint, 
the CLA that we are proposing is the first step of a procedure which eventually leads to the 
exact summation of the 3-body ladder series. Like perturbation theory and other approaches 
to the many-body problem, it allows systematic improvements, when required, at the cost of 
more computation. Physically, this simple approximation is well motivated and well balanced, 
and we provided evidence that it correctly describes the effect of the screening electron over 
the two final-state holes. The CLA is very accurate in the easy cases (high filling and/or small 
interaction); however, many transition metals are in the range of fillings and U where it is useful, 
and the approximation remains reasonable, Herglotz and qualitatively good even for U/W = 1. 
This holds true for both the one-body features and the two-hole-one-electron contributions 
which result from the full 3-body Green's function. The property of the CLA of remaining 
qualitatively correct even at rather strong coupling, when all other simple approaches badly 
fail, is an appreciable feature. It depends on the fact that the Core Approximation becomes 
accurate at strong coupling, which is the most critical regime for perturbation theory; on the 
other hand, it allows to carry on the theory to all orders, while treating the screening electron 
and the Auger holes on equal footing. We are currently working out the application of the new 
approach to the line shape analysis of actual experimental spectra. 

A The Auger and APECS currents 



In a simplified form of the One-Step model 20 of the CVV Auger spectra, the current of 



electrons with energy e& measured in a Auger-Electron-Spectroscopy (AES) experiment is 

poo /*oo 

J(e k )= dh dt2f{txM)e iek{tl ~ t2) - (A.l) 

Jo Jo 

Here, letting H$ represent the valence Hamiltonian without the core-hole, \ip) its ground state, 
H' the valence Hamiltonian in the presence of the core-hole, 

f(t u t 2 ) = ^(V'|e* ( ^^^|m)H ff i eiHs(tl ~ ta)fl 'iik / >W|e- i(fr '- ir) * 1 |V'), (A.2) 

m,m' 

where 

H A = y^ j M a 0a a a l3 . (A.3) 

a,/3 

produces the Auger holes with matrix elements M in spin-orbitals denoted by Greek symbols; 
r is an operator which produces virtual Auger transitions; core and free-electron operators have 
already been averaged out; the m, m' summations run over a complete set of valence states. 
In Ref. |Hj] a simple approximation was proposed. The basic idea was that the complete set 



of summations are largely exhausted by summing over just two orthogonal states, namely, 
\ip) and the relaxed initial state of the Auger transition, \<f>), that is, the ground state of the 
valence electrons in the presence of the core-hole potential. In this scheme, the Auger spectrum 

17 



has two main contributions, relaxed and unrelaxed. The unrelaxed contribution arises from 
m — m' — \ip), and can be expressed in terms of the two- hole Green's function G& . The 
relaxed contribution arises from m = m' = \<p), and is proportional to 

(<P\ala\e ms ^-^a pl a a ,\<P). (A.4) 

By a variational calculation Cini and Drchal plj showed that 



oc 

at 



J2 a l a ^)- ( A - 5 ) 



with a 1 creates an electron in a localized spin-orbital at the Auger site and ap annihilates an 
electron at the Fermi level. In this way the screening cloud is represented by a single electron 
that has moved from the Fermi surface to the empty local states of the Auger site. Using 
Equations ( |A.5|) and ( |A.4|) one obtains an expression for the relaxed contribution to the Auger 



line shape, involving the 3-body Green's function of Equation ( |1 . 1| ) . 

In APECS, on the other hand, one measures the Auger electron energy distribuition in 
concidence with the photoelectron energy, and one can selectively study the decay of the relaxed 



hole state. Semi-empirical studies [29| show that for a broad range of photoelectron kinetic 



energies the summation over m,m' is dominated by \(f>). 

B Proof of Equations (gHQ) and (g^TOf) 

In first quantisation the operator a^ap becomes \a)(j3\, and one can easily show that a^ (t) ap (£') 
becomes e lHot \a) (f3\e~ lHot . The average of any one-body operator on a Slater determinant is 
the sum of the averages on the occupied spin-orbitals k. Therefore, letting f k = 1 for occupied 
states and fk = for empty ones, 



(at 



(t)ap(t'))=Y,he^ {t - t ' ) (k\a)((3\k)- (B.1) 

k 

Therefore, summing over all the atomic spin-orbitals, one gets 

£>t (t)a, (t'))(a\(t')ap) = ^^f k f V e^^(k\a)^\k)e^^\ 7 ){l3\h^. (B.2) 

7 7 k,k' 

Now, exploiting the completeness of the 7 set, we obtain 

5>t ( t ) fl7 (*'))< a t it') ap ) = J2 f k e lekt (k\a) ((3\k>) = (oj, (t) a ). (B.3) 



In Equation ( B.3 ) the times t,t' are arbitrary. We are free to insert 6 functions to distinguish 



between positive and negative t and ensure that t' is intermediate between t and 0. This proves 
Equations ( p30|) . 
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C Application of the CLA 

From the diagonalisation of the one-body part of the problem, including the Hartree-Fock terms 
"ID, one readily obtains the propagator S ,e (ia,ja; t) where i — 1, 2 and a =f, |. Since there 



are two local orbitals in this model and the matrix {T^} is real the set of equations ( p. 16 ) 



represent a 3 x 3 linear problem that determines each component of the R^ 1 matrices, from 
which the B matrices Q2.25| ) can be easly obtained. 

For clarity, we report the system which determines R + : 

(-i) Sj (la, la; t) = R + (la, la, t) Sj (la, la; 0+) + R + (la, 2a, t) S^ (2a, la; 0+) , 
(-i) S£ (la, 2a; t) = R + (la, la, t) S% (la, 2a; 0+) + R+ (la, 2a, t) S* (2a, 2a; 0+) , (C.l) 
(-z) ^o (2o-, 2a; t) = R + (2a, la, t) S% (la, 2a; + ) + R + (2a, 2a, t) Sfi (2a, 2a; 0+) , 

yielding 

'R+ (la, la, t)\ fn n . a n X2 , a \ l ( S% (la, la, t)\ 

R + (la, 2a, t) = (-z) n 12 . a n 22 , CT S* (la, 2a, t) . (C.2) 

K R + (2a,2a,t)J \ n 12 , a n 22i J \S%(2a,2a,t)) 

where n^ )(r = Sq (ia,ja; + ). The diagonal n elements are much bigger than the off-diagonal, 
and the problem is well posed. A similar procedure gives the R~ functions using the propagators 
S^(ia,ja;t). 

We consider the following spin configurations: 

G (ijk, k'j'i'; u) = G (i t j | k T, k' T f 1 i' T; w) . (C.3) 

where the i,j,i',j' indices refer to the holes. To find them, we write down explicitly the CLA 
equations ( |2.24| ), where only the interactions involving the orbitals 1 and 2 at the Auger site 
s = 5 appear. In the Hubbard model the interaction terms are proportional to the product of 
number operators, and this entails the condition 

Upvpr = U8^ p 8 VjT . (C.4) 



That simplifies the set of equations ( |2.24| ) that reduce to a (8 x 8) linear system. It is convenient 



to use the shorthand notation where underlined numbers mark the position of R^ 1 functions, 
namely: 

Go (ijk, 111; w) = S* (i, 1; w) S% (j, 1; w) RT (k, 1; u) , 

G (ij_k, 111; w) = S* (z, 1; w) i? + (j, 1; w) S e (fc, 1; w) , (C.5) 

G (ijk, 111; w) = i2+ (i, 1; w) 5 h (j, 1; w) S e (A;, 1; w) , 
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then, we may write the system in the form 

G {ijk, k'j'i; u) - G CLA {ijk, k'j'i; u) = 

U{[G {ijk, 111; w) - Go {ijk, 111; w)] G (111, fc'/i'; w) + 
[G (ijfc, 211; w) - G (yife, 211; w) - G (£?&, 211; w)] G (112, fc'jV; w) + 

[Go (ul, 121; w) - G {ijk, 121; w)] G (121, fc'jV; w) + 
[G {ijk, 221; w) - G (yife, 221; w) - G (yifc, 221; w)] G (122, k'j'i'; u) + 
[Go {ijk, 112; w) - G (yife, 112; w) - G {ijk, 112; w)] G (211, k'j'i'; u) + 

[Go {ijk,212; to) - G (yifc, 212; w)] G (212, fc jY; w) + 
[Go {ijk, 122; to) - G (ylfe, 122; w) - G (yife, 122; w)] G (221, fc'jV; w) + 

[G {ijk, 222; w) - G (0A;, 222; w)] G (222, k'j'i'; u)} . (C.6) 

For the G sp contribution we use Equation ( |2.28|) for the propagator S {j [,j' [;uj), 

S {j i,j' Uu) = S {j i,j' i;u) -U 2 Y,Soti l,k l',w)V(k l,k' Uu) S {k' i,j' Uu) . (C.7) 

kk' 

Finally, we obtain the proper self-energies S {k l,k' 1;uj) from Equation ( j2.2£| ), where we insert 
the CLA form of the three body Green's function determined above. We recall that the one- 
body contributions arise from the anihilation of the up-spin electron and hole; therefore, one is 
left with the propagator for a down-spin hole, which is diagonal in spin because there are no 
spin- flip interactions in our model. Taking into account Equation flC.4|) , one finds 

s (l |, 1 1; uj) = G CLA (l T 1 1 1 T, l T 1 1 1 T; w) + G CLA (2 1 1 1 2 |, 2 1 1 1 2 |; u) + 

2G cla (1 T 1 1 1 | 7 2 t 1 I 2 T; cu) + G CLA (1 I 2 I 2 1, 2 I 2 1 1 j; w) , (C.8) 

S (1 |, 2 I; w) = G CLA (1 t 1 I 1 T, 1 T 2 I 1 T; w) + G CLA (1 | 1 I 1 T, 2 t 2 j 2 |; w) + 
G CLA (2 t U 2 T, 1 T 2 | 1 T; w) + G CLA (2 t 1 I 2 t, 2 T 2 I 2 T; cu) + 

G CLA (ll2j2l,ljll2I;u;), (C.9) 

£ (2 |, 2 |; to) = G CLA (1 t 2 | 1 t, 1 T 2 | 1 T; ^) + G CLA (2 | 2 | 2 |, 2 | 2 | 2 |; w) + 

2G cla (i | 2 | 1 t, 2 t 2 | 2 t; w) + G CLA (2 | 1 | 1 |, 1 | U 2 |; w) , (CIO) 



E(2i,U;w) = S(U,2|;w). (C.ll) 
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